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We develop a model of an ice shelf-ice stream system as a viscoelastic beam partially supported by an 
elastic foundation. When bedrock near the grounding line acts as a fulcrum, leverage from the ice shelf 
dropping at low tide can cause significant ( ~ 1 cm) uplift in the first few kilometers of grounded ice. 
This uplift and the corresponding depression at high tide lead to basal pressure variations of sufficient 
magnitude to influence subglacial hydrology. Tidal flexure may thus affect basal lubrication, sediment 
flow, and till strength, all of which are significant factors in ice-stream dynamics and grounding-line 
stability. Under certain circumstances, our results suggest the possibility of seawater being drawn into 
the subglacial water system. The presence of seawater beneath grounded ice would significantly 
change the radar reflectivity of the grounding zone and complicate the interpretation of grounded 
versus floating ice based on ice-penetrating radar observations. 

© 2012 Elsevier B.V. All rights reserved. 


1. Introduction 

Ice sheets transition from grounded streams to floating shelves 
over a several-kilometer-wide grounding zone that extends 
between the landward limit of tidally induced flexure and the 
hydrostatic line at which ice begins to oscillate with the full tidal 
range (Fig. 1). Grounding zones are critical to ice-sheet dynamics, 
because processes here determine the mass flux from grounded to 
floating ice, and thus eventually to the ocean. Furthermore, the 
interactions of ice-stream basal hydrological processes with the 
tidally changing ocean may exert important influences on water 
pressure, sliding and sediment flux near the grounding zone. 
Early studies focused on providing fundamental data on position 
of the grounding line (i.e„ the boundary between floating and 
grounded ice) and its motion in response to tidally induced 
flexure. Considerable success in this effort now motivates new 
studies examining the physical processes that take place at the 
ice-stream base in the grounding zone, and their effects on nearby 
subglacial hydrology. 

Flexure of ice shelves due to ocean tides has been studied 
at least since Robin (1958), who observed the phenomenon on 
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Maudheim Ice Shelf and noted that surface crevasses coincided 
with the location of maximum elastic stress. Holdsworth (1969) 
first modeled an ice shelf as an elastic beam, an approach 
followed in many subsequent studies (e.g., Holdsworth, 1977; 
Lingle et al„ 1981; Smith, 1991; Vaughan, 1995); Schmeltz et al. 
(2002) extended the model to a two-dimensional elastic plate. 
Several studies of this type could not reconcile observations with 
a single value of the elastic (Young’s) modulus for ice; Reeh et al. 
(2003) addressed this issue by adding viscoelastic behavior. 

These studies have been concerned mainly with the floating 
ice shelf, and often directed towards determining the location of 
the grounding line, a goal shared by many remote-sensing 
studies. Recent observations (e.g., Flicker et al., 2009; Brunt 
et al., 2010, 2011; Bindschadler et al., 2011; Rignot et al., 2011) 
have used radar and laser altimetry to precisely map grounding- 
line motion within grounding zones, building on earlier work by 
Scambos et al. (2007) with MODIS imagery. Sayag and Worster 
(2011) present interesting analytical and experimental results 
that help define the grounding-line position for an elastic ice 
shelf-ice sheet system at equilibrium, but do not yet include the 
effect of tides. 

Another main goal of geophysical studies, both observational 
and model-based, has been to understand the widely varying 
effects of tidal cycles on ice-stream flow (Bindschadler et al., 
2003; Anandakrishnan et al., 2003; Gudmundsson, 2007, 2011; 
Murray et al., 2007; Winberry et al., 2009, 2011; King et al., 2010; 
Walter et al., 201 1 ). The tides clearly are forcing strong variations 
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Fig. 1. Schematic of tidal flexure (not to scale). Inland limit of flexure denoted by F, 
grounding line by G, and hydrostatic line by H. Note that uplift inland of G is actually 
approximately 1% of the tidal height. 

in horizontal velocity as well as causing vertical viscoelastic 
bending, pointing to the need for a fully coupled solution. 

One important part of the path forward in developing such a 
solution comes from the tiltmeter studies of Smith (1991) and 
Jacobel and Bindschadler (1994), which show that significant 
bending can extend several km inland from the grounding line. 
Heinert and Riedel (2007) used a combination of GPS, gravi- 
meters, and tiltmeters to observe an upward displacement of 
1 cm approximately 1 km inland of the grounding line at Ekstroe- 
misen at low tide, implying the existence of a fulcrum near the 
grounding line. 

Studies addressing the ice-shelf flexural profile (e.g., Reeh et al., 
2003) have achieved considerable accuracy using models in which 
the tributary ice stream is “clamped" with zero deformation at and 
inland of the grounding line, because the tidal changes of order 1 m 
are large compared to the observed inland changes of order 1 cm 
(Heinert and Riedel, 2007; Jacobel and Bindschadler, 1994; Smith, 
1991). However, even the small inland uplift resulting from the 
leverage of the ice shelf as the tide falls could elastically transfer some 
of the weight of the ice, leading to a local drop in pressure beneath 
the ice stream. A sufficiently large pressure drop would affect ice- 
stream basal lubrication, till strength, and water exchange with the 
ocean, all of which are significant factors in grounding-line stability 
and ice-sheet mass balance. Thus, understanding the small tidal 
deflections inland of the grounding line seems to be a prerequisite 
for a fully coupled grounding-zone model. We therefore build on 
earlier work, but replace the clamp with a fulcrum at the grounding 
line (Fig. 1), to study this inland deformation. At this exploratory 
stage, we do not attempt the harder problem of allowing the 
grounding line to move. We note that the success of prior workers 
in quantitatively simulating the ice-shelf deformation assuming a 
clamped ice stream subject to the full tidal range indicates that 
grounding-line migration is small compared to the flexural wave- 
length of the ice shelf, likely because grounding lines tend to occupy 
relatively steep regions in basal topography (Anandakrishnan et al., 
2007), supporting the use of a fixed fulcrum. 

In the present study, we model tidal flexure of an ice shelf 
attached to an ice stream resting on an elastic foundation, with 
the goal of determining whether flexure-induced pressure varia- 
tions at the ice-stream bed are sufficiently large to influence 
subglacial water flow. 


2. Model derivation 

Our model is a viscoelastic modification of Euler-Bernoulli 
elastic beam bending (e.g., Turcotte and Schubert, 2002). When a 


thin beam (of any rheology) is bent, the balance of vertical forces 
on a small element leads to the relation 

8 x V=-q, (1) 

between the shear force V on a vertical cross section and the 
applied load q. The bending moment M is the integrated torque 
over a vertical cross-section of the beam: 


r0.5h 

M = ffxxZdz, (2) 

7-05/1 

where h is the beam thickness, a xx is the longitudinal stress, and 
the vertical coordinate z is zero on the beam centerline. The 
balance of torques then gives 

8 X M = V, (3) 


which can be substituted into (1) to obtain 

c$M = -q. (4) 


To find the bending equation for a material, its rheology must 
be substituted into the definition (2) of the bending moment. The 
Maxwell viscoelastic material is the simplest rheology that dis- 
plays both instantaneous elastic response and viscous damping. 
Reeh et al. (2003) argued that elastic-only models cannot capture 
all aspects of tidal flexure, while Gudmundsson (2011) demon- 
strated that the Maxwell model is sufficient to the purpose. 

For bending of a thin beam, the Maxwell rheology for the 
horizontal direction (the along-flow direction of the shelf) is given by 


£ 

1 -A 2 



(5) 


where £ is Young’s modulus, v is viscosity, e** is longitudinal strain, 
X is Poisson’s ratio, and dots indicate time derivatives. This may be 
derived from the equations for plane stress in a Maxwell material 
(Turcotte and Schubert, 2002) by assuming that the beam is thin in 
the vertical (z) direction and that bending occurs only in the xz-plane. 
In the differential operator form used by Reeh et al. (2003), and 
following his assumption that advective terms in the total time 
derivative are negligible, the Maxwell rheology may be written as 


Afrtxx) = 



2v(\-X 2 ) 


&XX — 



Exx = Q(6x*). 


( 6 ) 


Applying the operator P to the bending moment M results in a partial 
differential equation relating M and the vertical displacement w: 

dtM + ,2 M = -Ddtd^w, (7) 

2v(1 -/) 


where the flexural rigidity of the beam is defined as 


D = 


Eh 3 

12(1 -X 2 )' 


( 8 ) 


and we have used the curvature-based relation 

£** = —z8 x w, (9) 


between strain and vertical displacement (e.g., Turcotte and 
Schubert, 2002). 

Because we are interested in modeling flexure of the grounded 
ice as well as the ice shelf, we consider a partially supported 
beam. If the till underlying the ice stream is assumed to be elastic, 
the situation is analogous to a Winkler foundation in soil 
mechanics, and an extra term is added to (4): 

8^M = kw-q, (10) 

where the spring constant k for compression is Young’s modulus 
for the foundation divided by the thickness of the foundation 
(Jones, 1997). The applied load q is nonzero only for the floating 
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ice shelf, where it comes from the tide: 

Q = PswS( a sin cot—w), (11) 

where a and co are the tidal amplitude and frequency, and p sw is 
the density of seawater (Holdsworth, 1969). 

Finally, if we make the simplifying assumption that the 
viscosity v (but not the thickness h or flexural rigidity D) is 
constant with respect to x, we can combine the two previous 
equations for the bending moment M (taking c% of (7) and 
substituting (10)) to obtain 

Ek E 

d t (kw+c%(Dc%w))-\ ~-w = d t q-{ —q. (12) 

2v(l-2 2 ) 2v(l-2 2 ) 


We note that in the limit v->oo our viscoelastic bending 
equation reduces to the familiar equation for an elastic beam 
(e.g„ Turcotte and Schubert, 2002): 


dx 2 



q-kw. 


(13) 


The addition of viscous damping in our model affects the timing 
of the flexural response to the tide, while elastic behavior 
continues to dominate the solution. Our assumption of constant 
viscosity is thus reasonable for an exploratory model; c.f. the use 
of constant material parameters to successfully match observa- 
tions by Reeh et al. (2003). Model sensitivities to choices of 
parameter values are considered later in the paper. 

The equation is discretized and solved using Hermite cubic 
finite elements and fully implicit time differencing. Model resolu- 
tion is 100 m in space and 600 s in time. 


3. Simulations 

3.3. Control experiment 

The ice shelf-ice stream is 100 km long, with 50 km each of 
grounded and floating ice. The grounded ice has a uniform 
thickness of 1000 m, while the ice shelf tapers to 300 m thick at 
the ice front following a nonlinear profile similar to that used by 
Holland et al. (2008): 

h= — ti — » (14) 

(u 2 +(x-x gl )? /m 

where constants oq and a 2 are chosen to give the desired 
thicknesses at the grounding line and ice front, and x g/ is the 
position of the grounding line. We set m = 3 to obtain a profile 
closer to that of a shelf with negligible lateral drag (m=4) than to 
the linear profile (m = - 1) of a shelf in which lateral drag fully 
supports the driving stress. This relatively steep profile leads to a 
rapid decrease in rigidity (8), preventing the overestimate of 
leverage that would result from assuming an untapered shelf; 
some ice stream-ice shelf systems thus will have at least slightly 
more inland deflection than we calculate here. 

We assume that the upstream end of the ice stream is simply 
supported (w = 0, M = 0), while the vertical displacement at the 
ice front is given by the tide, and (consistent with observations) 
we impose 9 x w = 0 there. At the grounding line, we impose a 
fulcrum (w=0), so that the shelf acts as an overhanging beam 
(e.g., Hibbeler, 2000). We assume an elastic foundation under- 
lying the grounded ice, and test a wide range of spring constants 
to find those that are consistent with the ~ 1 cm maximum 
deflection observed by Heinert and Riedel (2007). We then check 
whether the range of consistent spring constants matches the 
range expected for possible basal conditions such as subglacial till 
over bedrock, and find good agreement. The situation we model 
is somewhat simpler than that of Heinert and Riedel (2007), 


who inferred transient grounding and ungrounding of Ekstromei- 
sen from a ridge downstream of the zero-tide grounding line, but 
for exploratory process modeling we do not believe that this 
should be a major issue. 

We set Young’s modulus and Poisson’s ratio for ice to 4.8 GPa 
and 0.4, respectively, following Gudmundsson (2011), who 
showed that these values make the behavior of a Maxwell 
material essentially identical to that of the more complicated 
Burgers material used by Reeh et al. (2003). We find that a 
viscosity of v=10 15 Pas allows for noticeable departure from 
purely elastic behavior while remaining high enough to prevent 
underdamping, consistent with observations. Model sensitivity to 
these and other parameters is discussed in the next section. 

Results for the control experiment, which uses a diurnal tide 
varying between + 1 m in height and a spring constant of 
k = 5 x 10 s Pa m -1 , are shown in Figs. 2 and 3. The grounded ice 
experiences a peak deflection of + 10.3 mm at the node 600 m 
inland of the grounding line, consistent with the 1 cm uplift 
reported by Heinert and Riedel (2007). Deflections greater than 



Fig. 2. Near-grounding line vertical deflection at 2-h intervals through one 24-h 
tidal cycle. Grounding line is at x=50 km. Rising tide (with positive height above 
datum) shown by x, falling (positive) by dashed line, falling (negative) by solid 
line, and rising (negative) by +. Each solution remains constant for x < 45 km or 
x > 70 km. 



Fig. 3. Closeup of Fig. 2 showing vertical deflection for grounded ice only. Note 
that the slight offset between corresponding deflections for rising and falling tides 
is caused by viscous damping. 
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1 mm continue up to 2.1 km inland of the grounding line, and 
deflections greater than 0.1 mm up to 4.3 km inland. 

On the ocean side, the majority of the shelf (at least the last 
30 km) remains very close to flotation, as expected. Deflections 
approximately 8 km offshore can exceed the tidal height by as 
much as 4.5 cm at high and low tide, a phenomenon observed in 
e.g., Reeh et al. (2003) which is similar to the formation of 
forebulges during flexure of lithospheric plates. At high tide, the 
deflection remains within 1 cm of the 1 m tidal amplitude starting 
12.0 km offshore of the grounding line, and within 1 mm starting 
17.3 km offshore. The slope of the deflection drops below 10~ 6 
(compared with a maximum of 2.15 x 10~ 4 ) starting 14.5 km 
offshore, and below 10~ 7 starting 20.2 km offshore. In observa- 
tional studies, the hydrostatic line is often estimated as “where 
the gradient of each elevation curve [the displacement w in our 
model] first tends to zero and is consistent with the tidal height 
predictions [our tidal forcing]” (Brunt et al., 2010); by this (partly 
visual) criterion, our hydrostatic line appears to be approximately 
14 km offshore. This distance is somewhat long relative to 
observations of the Ross and Filchner-Ronne Ice Shelves by 
Brunt et al. (2010, 2011); the difference is likely due to the 
steepness of our shelf profile, and possibly also to the high 
resolution of model output relative to satellite measurements. 

3.2. Sensitivities 

For each of the experiments described below, all parameters 
and model inputs other than those explicitly mentioned remain at 
their values for the control experiment. Experiments in which the 
domain length is extended (maintaining constant ice thickness in 
the extensions) do not produce any change in the near-grounding 
line solution, indicating that the results of the control experiment 
are not influenced by boundary effects. We focus here primarily 
on changes in the flexure of grounded ice resulting from variation 
of model parameters. 

Reducing Young’s modulus to 3.0 GPa and Poisson’s ratio to 
zero (thus reducing the rigidity D by over 47%) results in slightly 
greater maximum deflection ( + 10.4 mm) for the grounded ice, 
with a shortening of the area of nonzero deflection (maximum now 
only 500 m inland, and deflections greater than 1 mm only 1.7 km 
inland). Reducing the thickness of the grounded ice to 500 m 
decreases the rigidity by a factor of 8, leading to similar but larger 
changes (maximum deflection of +11.7 mm at 400 m inland, 
deflection greater than 1mm extending 1.2 km inland). These 
results are consistent with the general principle that response to 
an applied load becomes larger and more local as the rigidity of a 
beam decreases (e.g., Holdsworth, 1969). In contrast, changing to a 
linear shelf profile increases the rigidity of the lever while main- 
taining the rigidity of the grounded ice, resulting in the same area 
of nonzero deflection as the control experiment while increasing 
the maximum deflection to +11.8 mm. 

We find that reducing the viscosity significantly is not possible, 
as values below 3xl0 14 Pas cause the control experiment to 
become underdamped, with deflections growing over time. This 
underdamping is relatively weak for v ~ 10 14 Pa s, requiring several 
days for deflections to grow noticeably, but becomes stronger for 
v~10 13 Pas, with deflections greater than the elastic solution 
appearing within the first tidal cycle. Increased viscosity shortens 
the delay in response relative to a purely elastic beam, with the 
solution for v= 10 17 Pa s nearly indistinguishable from the elastic 
solution. We thus find that there is a relatively narrow range of 
values for which the viscosity significantly affects the model. 

As expected, the deflection (both offshore and inland) is a 
linear function of the tidal amplitude. The inland deflection is also 
strongly (but nonlinearly) dependent on the spring constant. 
When k is reduced by half, the maximum deflection increases 


to +14.2 mm at 700 m inland, with deflections greater than 
1 mm extending up to 2.6 km inland. Doubling k reduces the 
maximum deflection to + 7.4 mm at 500 m inland, with deflec- 
tions greater than 1 mm reaching only 1.6 km inland. The spring 
constant clearly plays an important role in limiting deflections to 
realistic values, as repeating the control experiment with k = 0 
(corresponding to the unphysical assumption that vertical motion 
of grounded ice has no effect on the pressure at its base) gives 
deflections greater than 2600 mm over 20 km inland. 

We now consider whether our empirically chosen value for the 
spring constant (k= 5 x 10 s Pa m -1 ) is physically plausible. We 
begin with the case of till overlying bedrock. With Young’s 
modulus for till ranging from 2 to 30 MPa (C. Marone, unpub- 
lished data from experiments in Rathbun et al. (2008); c.f. Dorf, 
2004), a typical layer thickness of a few meters leads to 
10 7 Pam~' as a reasonable value for k in the till. Similarly, 
bedrock with Young’s modulus of 0(10) GPa and elastic response 
likely extending a few km deep (comparable to the flexural 
wavelength of the ice) leads to fc= 10 7 Pa lrr 1 . The combined 
elastic response may be modeled as a compound spring consisting 
of two springs connected in series, so that each spring bears the 
full load. The spring constant is then given by k = krk 2 /(k\ + k 2 ) 
(e.g., Nautiyal, 2001), which in this case gives our chosen value of 
k = 5 x 10 6 Pa m 1 . Considering reasonable values for both mate- 
rial properties and observed deflections suggests that k should lie 
roughly in the range of 1 x 10 6 to 3 x 10 7 Pa m -1 , corresponding 
to maximum deflections between 0.44 and 2.2 cm for our control 
experiment. We also consider the addition of a thin layer (a few 
cm) of water, assuming for simplicity that the till is sufficiently 
saturated for this case to be modeled as a three-component 
compound spring. As the bulk modulus of water is over 2 GPa, 
the spring constant for this layer will be on the order of 
10 n Pa m _1 . This is so much stiffer than either the till or bedrock 
layer that its addition in series would have (for our purposes) 
negligible effect on the spring constant of the till-bedrock 
combination. 


4. Implications for subglacial hydrology 

Because we have an elastic foundation underlying the 
grounded ice, the change in subglacial pressure at any point on 
the ice stream is just AP= -kw. (We thus implicitly assume that 
water, where present, is locally redistributed quickly enough to 
avoid formation of air-filled cavities; however, this does not 
yet allow extensive horizontal flow of water in response to the 
changing pressure.) Applying this to our control experiment, we 
see that the maximum pressure drop beneath grounded ice 
induced by low tide is 51.5 kPa, located 600 m inland at the point 
of maximum uplift; we note that an average density of 
1028 kg m -3 for seawater gives a pressure drop of only 
10.08 kPa in the ocean when the tide falls by 1 m. The pressure 
gradient beneath grounded ice for the control run at low tide, 
from the calculated distribution of uplift with spatially constant k 
(Fig. 4), ranges from 200 Pa m -1 at the grounding line to 
42 Pa m _1 in the opposite direction at 1.2 km inland. 

To assess the significance of this pressure gradient, we com- 
pare it with typical (zero-flexure) values of the gradient of the 
subglacial hydraulic potential (f> h : 

V<t> h =PigVS+(p w - Pi )gS/B, (15) 

where S and B are the surface and basal elevations of the ice 
stream (Shreve, 1972; Cuffey and Paterson, 2010). For ice streams 
similar to those on the Siple Coast, with surface slopes of 5 x 10~ 3 
or less, V</> h is generally below 50Pam _1 ; on near-grounding- 
line ice plains, the hydraulic gradient can be much lower. We thus 
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Fig. 4. Pressure gradient under grounded ice for control experiment at high (solid 
line) and low (dashed line) tide. 


conclude that the pressure gradient caused by ice pressure could 
significantly influence subglacial water flow, although further 
investigation of the detailed structure of the water system is 
necessary to fully assess its effect. In particular, the relative sizes 
of the flexural pressure gradient and hydraulic gradient suggest 
the possibility of ocean water being drawn into the subglacial 
water system near low tide, with a brackish mixture being forced 
farther inland when the pressure gradient reverses near high tide. 

The maximum pressure drop generated in the model is clearly 
dependent on the spring constant k. To investigate this sensitivity, 
we calculate the maximum pressure drop and corresponding 
uplift at low tide for a broad range of k (Fig. 5). The nonlinear 
relationship between k and max(w) leads to greater uplift for 
smaller pressure drops, which is consistent with interpreting the 
pressure drop as the self-weight that must be supported by the 
ice. For the model configuration used in the control experiment, 
the subglacial pressure drop exceeds the ocean pressure drop 
(10.08 kPa, as noted above) when the maximum uplift is less than 
4 cm; this difference becomes large as the maximum uplift falls 
toward the observed value of approximately 1 cm. This suggests 
that ocean water being drawn beneath the ice stream is possible 
for a broad range of k, again pointing to a need for further 
investigation. 


5. Discussion 

Observational confirmation of inland flexure is difficult due to 
the low amplitude of the vertical displacement. Vertical changes 
of 1 cm cannot be reliably measured using satellite laser altimetry 
(Flicker et al., 2005) and are at the resolution limits of GPS on 
floating ice (King, 2004). Differential interferometric synthetic 
aperture radar is also unable to directly detect the inland flexure 
limit F (Fig. 1); studies relying on this method (Rignot, 1998; 
Rignot et al., 201 1 ) have calculated F by fitting an elastic-beam 
model to observations. Other geophysical techniques, however, 
namely tiltmeters for surface observations and ice-penetrating 
radar (IPR) for basal observations, show more promise. Tiltmeters 
generally have sub-millimeter accuracy (Smith, 1991) and have 
been used to locate the limit of inland flexure along grounding 
lines in the Ross and Ronne-Filchner embayments (Smith, 1991; 
Jacobel and Bindschadler, 1994). Recent ground-based ice-pene- 
trating radar observations along the grounding line of West 



Fig. 5. Pressure drop as a function of maximum uplift of grounded ice at low tide. 
Note that 10.08 kPa corresponds to the pressure drop in the ocean at a tidal height 
of - 1 m. Except for the spring constant k, model parameters are as for the control 
experiment. 


Antarctica’s Siple Coast and East Antarctica’s Dronning Maud 
Land, which rigorously account for spatial variation in englacial 
attenuation using direct measurements (MacGregor et al., 2011) 
or multiple models (Matsuoka et al., 2012), indicate that ground- 
ing lines, as depicted on ice-penetrating radar, are spatially 
variable. On the Siple Coast, lack of a change in amplitude of 
the basal reflector across the grounding line indicates that there 
may be a brackish layer of water well upstream ( ~ 10 km of the 
grounding line) (MacGregor et al., 2011); in contrast, observations 
in Dronning Maud Land show a sharp contrast in basal reflectiv- 
ity, indicating the grounding line here may be a barrier to 
seawater intrusion (Matsuoka et al„ 2012). A complete investiga- 
tion of large, existing airborne IPR datasets could indicate areas 
where seawater intrusion may occur inland of grounding lines 
and whether these areas are correlated with areas where tidal 
flexure may locally affect hydropotential enough to force sea- 
water infiltration upstream of grounding lines. Access to the 
upstream and downstream parts of the grounding zone, by bore- 
hole and submersible, would be of great interest and may be 
required to gain full understanding. 

The possibility of seawater intrusion strongly suggests that 
future work on near-grounding-line tidal flexure will require 
coupling a model of subglacial hydrology with a viscoelastic 
beam model. Even for the fairly large pressure gradients induced 
by ice flexure, it is unlikely that Darcian flow through the till 
(e.g., Murray and Clarke, 1995; Pimentel and Flowers, 2011) could 
transport enough water during a partial tidal cycle to greatly 
reduce the pressure gradient, given reasonable values of effective 
conductivity and layer thickness. This suggests that the flexure 
will affect hydrologic potential and thus flow in the subglacial 
water system, whether it consists of a sheet (Weertman, 1972), 
a network of shallow canals (Walder and Fowler, 1994), or 
cavities and channels in some combination (Schoof, 2010). In 
turn, this is likely to affect lubrication of ice flow, and sediment 
transport. Modeling this system will be complicated by the 
requirement that the water discharge to the ocean at the onset 
of flotation, rather than flowing at a water pressure slightly less 
than the ice overburden and controlled by local ice-water- 
sediment interactions, as occurs well upstream of the grounding 
zone. Further complexity will be introduced by the effects 
of changing load and water pressure on till properties, and 
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additional modeling may need to consider a fulcrum that moves 
during the tidal cycle. We anticipate that considerable modeling 
and observational work will be required to understand this 
system fully and model it accurately across the range of possible 
conditions. 

If a significant fraction of the peak uplift of grounded ice at low 
tide is accommodated by lateral inflow of water, then both the 
pressure drop and the elastic deformation of the substrate there 
will be smaller than calculated above; however, the pressure drop 
must remain significant to drive the water inflow. Given the 
known wide range of hydrological systems beneath glacial ice 
(e.g, Walder and Fowler, 1994; Schoof, 2010; Pimentel and 
Flowers, 2011), diverse behavior at different grounding lines 
seems likely. 

The tidal flexure of the downglacier part of the ice stream will 
load and unload any till beneath. A range of behaviors is possible, 
but with a general tendency for cyclic loading at such frequencies 
to favor net compaction (e.g., Murray and Clarke, 1995; Lekarp 
et al., 2000). If this proves to apply to an ice stream that is 
lubricated by a soft, deforming till, then the glacier bed in the 
region so affected will be relatively stiff, slowing the ice flow and 
stabilizing the grounding line. This is consistent with observa- 
tions showing that the surface of the Siple Coast ice streams steps 
down relatively steeply in crossing a relatively stable grounding- 
line position into the Ross Ice Shelf (Florgan and Anandakrishnan, 
2006), and with the GPS and seismic evidence for a “sticky spot" 
acting near the grounding zone of one of these ice streams 
(Whillans) during its twice-daily events of rapid slip (Winberry 
et al„ 2011). These are clear targets for additional field studies. 

6. Conclusions 

Physical understanding and the limited available data indicate 
that tidal flexure of ice-stream/ice-shelf systems should cause 
vertical displacements extending inland of the grounding line. Our 
viscoelastic flowline model, including elastic response of likely 
subglacial materials, produces vertical displacement of order 1 cm 
on the order of 1 km inland, consistent with available data. The 
associated subglacial pressure variations are larger than the tidal 
effect in the free ocean, with gradients along flow that reach peak 
values much larger than those farther inland driving subglacial 
water flow toward the grounding zone. We anticipate that addi- 
tional modeling of the fully coupled system will show that the large, 
oscillating subglacial pressure gradients in and near the grounding 
zone exert important controls on the coupled flow of ice, water and 
till, perhaps including transport of seawater well inland. Further 
observations with ice-penetrating radar will be needed to locate 
areas of seawater intrusion beneath grounded ice and assess their 
effects on grounding-line detection by radar reflectivity. 
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